Non-Adiabatic Fluctuation in Measured Geometric Phase 
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We study how the non-adiabatic effect causes the observabfe fluctuation in the "geometric phase" 
for a two-levei system, which is defined as the experimentaily measurabie quantity in the adiabatic 
limit. From the Rabi's exact solution to this model, we give a reasonable explanation to the ex- 
perimental discovery of phase fluctuation in the superconducting circuit system [P. J. Leek, et ai, 
Science 318, 1889 (2007)], which seemed to be regarded as the conventional experimental error. 
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Introduction - It was discovered by I. I. Rabi [l[ that 
the non-adiabatic transition of a quantum system in a 
time-dependent magnetic field was subject to the sign 
of its magnetic momentum. The exact solution was first 
given in 1937, but its physical significance for the relative 
phase acquired under adiabatic evolution was not clari- 
fied until five decades later 0. It was M. V. Berry who 
found that this phase might contain a geometric part, 
now called the Berry's phase. Then the quantum adia- 
batic approximation theorem (QAAT) Of was reproved 
to naturally include the Berry's phase [4[ and general- 
ized to deal with the non-adiabatic effects for many cases 
[EH, 0,11) Hi- On the other hand, because of its geometric 
dependence, conditional geometric phase was proposed as 
an intrinsically fault-tolerant way of performing quantum 
computation [111 ]. 

In this paper, associated with a recent experiment 
about phase fluctuation in the superconducting circuit 
system [3], the above Rabi's solution is used to study 
in details the non-adiabatic effects for a two-level system 
(TLS) in a harmonically rotated field (see Fig.l). This 
field can be realized with a microwave field perpendicular 
to the static magnetic field both applied to the system. 
With the phase of the microwave linearly varying with 
time, the Hamiltonian harmonically rotates in the para- 
metric space. 

Generally speaking, the Berry's phase is always accom- 
panied with the dynamical phase, and thus its pure ef- 
fect can not be observed directly. However, we apply a 
7r-pulse to the TLS so that the evolution is divided into 
two parts with both of them in the same path but in the 
opposite directions. In this case, the effect of the dy- 
namical phase can be completely eliminated. This is the 
technique referred to spin echo technique [l3[ . Thus, the 
pure geometric effect can be observed. And that may re- 
sult in observable fluctuation of the measured geometric 
phase due to the Rabi's non-adiabatic transitions. 

Non-Adiabatic Effect with Berry's Phase.-The evolu- 
tion of the system can be well described with the Hamil- 
tonian 

H(t) = ^(Ao- z + tt R a x cos uj R t + Q R a y sin iv R t), (1) 
where A is the energy splitting without the microwave 



static magnetic field 





microwave field 



FIG. 1: (color online) (a) Schematic diagram of a TLS in a 
static magnetic field along with a microwave field, (b) Real- 
ized effective Hamiltonian rotating in the parametric space. 



field, rt R the Rabi frequency of the microwave, tu R the 
oscillating frequency of the microwave's phase, o~ x>v , z the 
Pauli matrixes. Note that Hamiltonian QT) is exactly 
the effective Hamiltonian realized in Ref.|l4|, where the 
rotating wave approximation [l5| was applied. Straight- 
forwardly, its instantaneous eigen states are obtained as 



|e(f)>=cos-|0)+sm-e^|l), (2) 
\g(t))=sm 6 -e- i » «*|0) - cob ||1>, (3) 



Here the en- 
and the mixing angle 



with corresponding eigen energies ±w/2 
ergy splitting is ui — -\/A 2 + Q R , 

is 6 = tan -1 (S7j,;/A). We also emphasize that due to the 
requirement of the single valueness of the eigenfunctions 
for a given Hamiltonian without singularity, the phase 
factor exp(±iu} R t) in \e(t)} or \g(t)) is fixed once the fac- 
tor in the other state is chosen [L2I ]. 

At time t, the evolution state is assumed to be a su- 
perposition \tp(t)) — a(t)\e(t)} + (3{t)\g(t)) of two instan- 
taneous eigen states . The time-dependent Schrodingcr 
equation H\^f(i)) = idt \^{t)) leads to the following equa- 
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tions of coefficients 

a = + ujB.sm 2 -)a + sin 6, (4) 

P' = - lur cos 2 -)/3' + ia^-sin9, (5) 

where (3'{t) — fi{t) exp(-iojRt). 
Under the adiabatic conditions 

| + WflS in 2 (-)»^sm0, (6) 
- - ujr cos (-) > — sm 9, (7) 

the adiabatic approximate solutions to Eqs.(|4|) and (O is 
obtained by ignoring the terms with LURsin9 /2 . They 
show that both norms of the amplitudes remain the same 
as their initial values, while they acquire a Berry's geo- 
metric phase ±u>Rt(l — cos#)/2 in addition to the dy- 
namical phase ±wt/2 respectively. 

On the other hand, the above Eqs.([3]) and can be 
solved exactly and it should be done so when the adi- 
abatic condition is broken under certain circumstances. 
Thus, we have 

a = A 1 e iu, + t +A 2 e iw - t , (8) 
= B ie luJ+t + B 2 e lu> -\ (9) 

where u>± = (— u>r ± f2)/2 = (— ljr ± 
\/ ijj 2 — 2<jjujr cos 9 + u> R )/2, and the coefficients are 
determined by the initial values a(0) and /3(0), 

_ P(0)tURsm9 + a(0)(-uj + ^ + ) 
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sin(9 + a(0)(u; + £_) 
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S_) + a(0)ujR sin 9 
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+ S+) - a(0)wR sin6» 


2Vt 



where E± = lur{1 ± cos0) + 2lo + . 

Generally speaking, the Berry's phase can not be ob- 
served directly from the experiment, since the dynamical 
phase always occurs along with the Berry's phase. Here, 
we consider a concrete case where the total evolution is 
divided into two rounds, both of which are in the same 
path but with opposite directions. Provided that both 
the excited and the ground states acquire the same dy- 
namical phase but with opposite signs in each round, the 
dynamical phase can be canceled when the amplitudes 
are exchanged after the first round of evolution. There- 
after, by experiencing the inverse rotation in the para- 
metric space, the dynamical phase is canceled while the 
Berry's phase is doubled since the latter depends on the 
sign of the angle velocity uir while the former does not. 
Excellent agreement may be expected with the Berry's 
prediction 

4>b = 2uj r T{\- cos 9) (10) 
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FIG. 2: (color online) (a) Comparison between the Berry's 
phase (blue solid line) (j>B and non-adiabatic phase (red cross) 
4>na with n — 1 circular rotation in each round, (b) The dis- 
crepancy between them with red solid line for the numerical 
result and black dashed line for the second order approxima- 
tion. 

provided that the adiabatic condition, i.e., Eq.©, is sat- 
isfied. Here, T is the evolution time for each round. 

However, when the non-adiabatic effect is considered, 
small deviation is expected. In Fig. ([2]), we plot the exact 
phase calculated from Eqs.© and ©, which is defined 
as the phase of a(2T)/3*(2T), denoted as 

Ka = angle(a(2T)/r(2T)). (11) 

This is an observable quantity in experiment, which can 
be determined by measuring the complex amplitudes 
a(2T) and (3(2T). It just recovers the Berry's phase in 
adiabatic limit. 

It is predicted that the Berry's phase is proportional 
to the solid angle = 2ir(l — cos 9) subtended by the 
path. Recently, a measurement of the Berry's phase in 
the superconducting qubit was carried out [14| . In order 
to compare the theoretical analysis with the experimen- 
tal result, we adopt the same parameters as those given 
in Ref.pJ], i-e., A = 50 MHz, oj r = (4n+l) MHz with n 
being the number of loops. The tiny difference between 
those two can almost not be distinguished in Fig. ([2k). 
Moreover, as shown in Fig.([2jD), there are small oscilla- 
tions in the deviation between them, Acf> = <p na — 4>b, 
with the root-mean-square deviation of 0.015 rad from 
the expected lines while the counterpart for n — 1.5 is 
0.043 rad. They are in reasonable agreement with the ex- 
perimental result, i.e., 0.14 rad [141 . considering that the 
rotating wave approximation [l5| was applied to obtain 
the Hamiltonian (fTJ and the non-adiabatic effect in the 
process of applying the microwave field also accounted 
for part of the deviation. 

Second-Order Fluctuation. -It was C. N. Yang who 
firstly pointed out that the Berry's phase could be re- 
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covered from the original QAAT by retaining the first 
order term 0(lo r /lu) in the phase. And this point of 
view was shortly confirmed by one (CPS) of the authors 
@ using the Rabi's exact solution. Here, with careful 
calculation to the second order term 0(uj r /oj) 2 , the fluc- 
tuation in the phase is obtained. For an initial state with 
a(0) = /3(0) = l/v2) the measured deviation from the 
expected Berry's phase is given as 

A0 = A 2 (sin 0i + 4 sin 02 + 2 sin 03 + 2 sin 04 
+ 2 sin 05 + sin 06 + 2 sin 07 + sin 0g 
+ 2sin0 9 + 2sin0io + 4sin0n + 4sin0i 2 ), (12) 

where A = lo r sin 9/2lu, <j)j = 0^- — 4> R (J = 1,2 . . .) with 

01 = 7r- 2Tlu r cos9, 

02 = -TUJR, 

03 =Tw R (-l + 2cos0), 

, T(uj 2 r + 8lurlu + Acu 2 — Auj r uj cos 9 — uj 2 r cos 29) 



t>12 



Aui 

T(-(vr + 2uj) 2 + 4u R u cos + ui R cos 20) 

" + 4^ ' 

T((wr-2w) 2 - wfjcos20) 
2lu 

T(uj 2 r - 2u! R uj + Au? - uj 2 r cos 20) 
2^ ' 
T{lj 2 r + Auj 2 -lj 2 r cos 20) 

" + 2^ ' 

T((w fl - 2oj) 2 - Au r lo cos 9 - wfj cos 20) 
7F+ 4^ ' 
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It can be seen from Eq. (fT2|) that the sinusoidal devi- 
ation from the expectation will be present in the mea- 
sured phase. As shown in FigJ2Jb), the behavior is well 
described by the second order approximation. Although 
there is small distinction between it and the exact result, 
it is believed that it is due to the higher order terms. 
Based on the above theoretical analysis, we may safely 
arrive at the conclusion that due to the non-adiabatic 
evolution of the TLS, not only do the relative phases of 
the amplitudes of the two states change with time, but 
also do their norms. In other words, the non-adiabatic 
transition results in the small fluctuation of the measured 
phase. 

To further explore the non-adiabatic effect, we investi- 
gate the phase fluctuation for different rotation velocity 
u)rS. In Fig|3l the discrepancy between the measured 
phase and the Berry's phase is plotted. As expected from 
Eq. (fT!?)) . the sinusoidal oscillation is again witnessed. Ad- 
ditionally, the oscillating amplitudes rise as the rotating 
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FIG. 3: (color online) Comparison between different numbers 
of rotation loops n, with black dotted line for n = 1, blue 
dashed line for n = 2, red dash-dot line for n — 3 and green 
solid line for n — 4. 



velocity is increased. Here, the total time for evolution 
is fixed [l(| • This result is consistent with Eq. (fi"2|) as A0 
scales as u>\. Notice that all of them nonexceptionally 
approach zero at the both ends. It is a reasonable result 
since A0 vanishes as = 0. Actually, the Hamiltonian 
remains the same as its initial state in the parametric 
space. The only effect for time evolution is to acquire a 
dynamical phase. On the other hand, as approaches 
tt/2, 



uj = J tt 2 R + A 2 > fl R > A > lo r 

since A is fixed. Due to vanishing of A, we have a zero 
deviation from the predicted phase, i.e., A0 = 0. In the 
limit -> tt/2, the initial state |V>(0)) = (|0) + |1))/V2 
is the eigen state of the initial Hamiltonian H(0) ~ 
Q R a x /2. Intuitionally, the evolution of the system is sim- 
ilar to the situation that a classical magnetic moment 
initially parallel to the applied field closely follows the 
rotation of the field. 

Conclusion. -hi this paper, we have investigated the 
fluctuation in the phase due to the non-adiabatic evo- 
lution. In contrast to the adiabatic evolution, the si- 
nusoidal deviation from the expected line drawn for the 
Berry's phase is observed. The phase fluctuations discov- 
ered in the recent experiment can be partially explained 
by our theoretical analysis. For a given time of evolu- 
tion, it is predicted that the fluctuation from the Berry's 
prediction rises larger and larger as the number of loops 
is increased. 
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